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ABSTRACT 

We investigate standing kink magnetohydrodynamic (MHD) oscillations in a 
prominence fine structure modeled as a straight and cylindrical magnetic tube only 
partially filled with the prominence material, and with its ends fixed at two rigid walls 



P_i. representing the solar photosphere. The prominence plasma is partially ionized and 

Q ! a transverse inhomogeneous transitional layer is included between the prominence 

+i ! thread and the coronal medium. Thus, ion-neutral collisions and resonant absorption 

c^ . are the considered damping mechanisms. Approximate analytical expressions of the 

period, the damping time, and their ratio are derived for the fundamental mode in 

^ ! the thin tube and thin boundary approximations. We find that the dominant damp- 

0> ! ing mechanism is resonant absorption, which provides damping ratios in agreement 

Q>^ ! with the observations, whereas ion-neutral collisions are irrelevant for the damping. 

The values of the damping ratio are independent of both the prominence thread length 

and its position within the magnetic tube, and coincide with the values for a tube fully 

O ! filled with the prominence plasma. The implications of our results in the context of the 

MHD seismology technique are discussed, pointing out that the reported short-period 

(2 - 10 min) and short-wavelength (700 - 8,000 km) thread oscillations may not be 

rS [ consistent with a standing mode interpretation and could be related to propagating 

5^ \ waves. Finally, we show that the inversion of some prominence physical parameters, 

e.g., Alfven speed, magnetic field strength, transverse inhomogeneity length-scale, 

etc., is possible using observationally determined values of the period and damping 

time of the oscillations along with the analytical approximations of these quantities. 



Subject headings: Sun: oscillations — Sun: corona — Sun: magnetic fields — waves 
— prominences/filaments 
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1. INTRODUCTION 



Oscillations and propagating waves are com monly reported in observations of solar promi 



nences and filaments (see recent reviews by, e.g., lBallestedl2006l : lEngvoldl l2008l : iMackay et al. 



2010h . In high-resolution observations, transverse oscillations of prominence fine structures are 
frequently detected. These fine structures, here called threads, appear as a myriad of long (5" -20") 
and thin (0" . 2 - ".6) dark ribbons in Ha images of filaments on the solar disk (e.g., iLin et al. 



2007L I2OO8L |2009|) . as well as in observations of prominenc es in the solar limb from the So 



lar Optical Telescope ( S OT) aboard the H inode satellite (e.g.. lOkamoto et al.l 120071 : iBerger et al 



20081 : IChae et al.l l2008l : iNing et al.ll2009X From the theoretical point of view, prominence fine 



structures have been rnodeled as magnetic flux tubes anchored in the solar photosphere (e.g., 
Ballester & Priestlll989l : iRempel et al.lll999h . which are piled up to form the prominence body. 
In this interpretation, only part of the flux tubes would be filled with the cool (~ 10^ K) filament 
material, which would correspond to the observed threads, while the rest of the magnetic tube, i.e., 
the so-called evacuated zone, would be occupied by hot coronal plasma. 

Common features of the transverse oscillations of prominence fine structures detected in 
Doppler signals and Ha sequences are that the reported periods are usually in a narrow range 
between 2 and 10 minutes, that the velocity amplitudes are smaller than ~ 3 km s'\ and that the 
oscillations seem to be damped after a few periods. Typically, the n umber of oscillatory periods 
observed before the oscillations disappear is les s than 10 (see, e.g.. iMolowny-Horas et al.lll999l : 
Terradas et al.ll2002l : lLinll2004l : iNing et al.l 120091) . Theoretically, the oscillations have been inter- 
preted in terms of kink magnetohydrodynamic (MHD) modes supp orted by the fine stru cture, mod 



eled as a magnetic slab (iJoarder et al.lll997 



2002Ll2010l : lDvmova & RudermanI 



2005 



Diaz et al.ll2001 



Terradas et al. 



2003!) or a cylindrical tube (IDiaz et al. 



20081) partially filled with the prominence 



plasma, whereas several damping mechanisms have been prop osed to explain the quick attenuation 
(see, e.g.. loiiveijl2009l : lArregui & Ballesteill201ol : ISoleiJl2oTol) . 



By neglecting the variation of the plasma parameters along the fine structure and adopting a 
prominence thread model composed of a ho mogeneous rnagnet ic flux tube with prominence con- 
ditions embedded in a coronal environment, ISoler et al.l (|2008|) studied the temporal damping of 
propagating kink MHD w aves due to nonadia batic eff'ects (radiative losses, thermal conduction, 
and plasma heating), while ISoler et al.l (l2009b|) investigated the attenuation in the same configura- 
tion but considering ion-neutral collisions as the damping mechanism. These authors concluded 
that neither nonadiabatic effects nor ion-neutral collis ions can produce kink mode damping times 
compatible with those observed. On the other hand, lArregui et al.l (l2008b|) considered a similar 
model but neglected gas pressure (i.e., the yS = approximation, with /3 the ratio of the gas pres- 
sure to the magnetic pressure) and took into account the presence of a transversely inhomogeneous 
transitional layer between the thread and the coronal plasma. In such a configuration, the kink 
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mode is resonantly coupled to Alfven continuum mo des, and so the kink mode is damped by the 
process of resonant absorption. lArregui et alj (|2008bh numerically obtained values of the damping 
time that are consistent with those reported in the observations. Resonant absorption has been 
previously propos ed as an explanation for both the temporal dampin g of coronal loop transverse 
oscillations (e.g.. iRuderman & RobertsI 2002 : ^Goos sens et al. 12002') and the spatial damping of 



propaga ting kink waves (e.g., 



Pascoe et al. 



2010.:.Terradasetal. 



2010 ). Subsequently. ISoler et al. 



(l2009al ) performed a more in-depth analytical and numerical investigation of the damping by res- 
onant ab sorption in prominenc e threads by inclu ding gas pressure, a nd obtained similar results to 
those of lArregui et al.l (|2008bf) . More recently, ISoler et al.l (|2009cl) studied the combined effect 
of resonant absorption and the promin ence plasma partial ionization on the kink mode damping. 
For realistic values of the wavelength, ISoler et al.l (|2009c|) concluded that partial ionization does 
not affect the process of resonant absorption, and so the obtained value s of the daniping ti me are 
the same as in a fully ionized prominence thread (^Arre gui et al.ll2008bl : ISoler et al.ll2009a|) . Ion- 
neutral collisions may become more relevant than resonant absorption for the temporal damping 
of propagating kink modes when short wavelengths of the order of 10^ km and smaller values are 
considere d, while the range of typi cally observed wavelengths in prominences is between 5x10^ 
- 10^ km doiiver & Ballesteill2002l). 



Whether the reported observations of transversely oscillating filament and prominence threads 
are related to p ropagating waves or standing oscillations is a subject of debate. For example, 
Lin et al.l(|2009l) explained their observations in terms of propagating kink waves, w hereas iTerradaset al 
2008h proposed standing oscillations as an explanation of the observations by lOkamoto et al 



(|2007|) . It is likely that propagating waves may be generated by localized disturbances in, e.g., the 
footpoints of the magnetic tube, while standing oscillations may be related to more global pertur- 
bations of the whole magnetic structure. Regarding the damping, all the works cited above studied 
the temporal damping of propagating kink waves, hence the wavelength (or the wavenumber) is a 
free parameter in their case, but the problem of the damping of standing oscillations has not been 
addressed yet in the context of prominence fine structures. Moreover, the effect on the damping 
of the longitudinal variation of the plasma parameters along the fine structure was not addressed 
in these previous works. Therefore, the aim of the present investigation is to broach the problem 
of the damping of standing kink MHD oscillations of longitudinally nonuniform prominence fine 
structures. 



The model config uration adopt ed here is simi l ar to th at considered by lDiaz et al.l(|2002ll2010|) . 
Dymova & RudermanI (|2005|) . and iTerradas et al.1 (l2008b . namely a straight and cylindrical mag- 
netic flux tube only partially filled with prominence plasma. The rest of the magnetic tube, as 
well as the external medium, has typical coronal properties. The/? = approximation is adopted 
for the sake of simplicity since we restrict ourselves to kink modes, which are correctly described 
by this approximation. Standing oscillations are studied by imposing the line-tying condition at 



the ends of the cylinder. As in ISoler et alj (|2009c|) . we assume that the prominence plasma is 
partially ionized and include a transversely inhomogeneous transitional layer between the dense 
prominence thread and the external corona. Hence, the mechanisms of ion-neutral collisions and 
res onant absorption are conside red as damping mechanisms. We follow the method introduced 
by iDymova & RudermanI (120051) based on the thin tube limit, and derive a dispersion relation for 
damped kink MHD oscillations in the thin boundary approximation. Analytical expressions of 
the period, the damping time, and the ratio of the damping time to the period are obtained, while 
a general parametri c study is performed by numerically solving the full dispersion relation. In 
a subsequent work, lArregui et al.l (120101) investigate the damping of kink oscillations beyond the 
thin tube and thin boundary approximations by numerically solving the full resistive eigenvalue 
problem in two-dimensional, nonuniform threads. 

This paper is organized as follows. Section [2] includes a description of the model configura- 
tion and the mathematical method. The dispersion relation of standing kink MHD oscillations is 
obtained in Section [31 which also contains analytical approximations. The results of solving the 
dispersion relation are given in Section HI while their implications for prominence seismology are 
discussed in Section [5] Finally, Section [6] contains our main conclusions. 



2. MODEL AND METHOD 



2.1. Equilibrium Configuration 
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Fig. 1 . — Sketch of the model configuration adopted in this work. 



The model considered here is schematically plotted in Figure \T\ We consider a straight and 
cylindrical magnetic tube of length L and radius a, whose ends are fixed by two rigid walls repre- 



senting the solar photosphere. The magnetic tube is only partially filled with the cool and dense 
prominence material, and is composed of a dense region of length Lp with prominence conditions 
and representing the prominence thread, surrounded by two much less dense zones corresponding 
to the evacuated part of the tube. According to the observed typical values o f thread widths and 
lengths from the high-resolution observations (e.g.. lLinll2004l : lLin et al.ll2008|) . the ranges of real- 
istic values of a and Lp are 50 km < a < 300 km and 3,000 km < Lp < 28,000 km. On the other 
hand, the total tube length, L, cannot be measured from the observations, but one can relate L to 
the typical spatial scale in prominences and filaments, i.e., L ~ 10^ km. 

For simplicity, both the prominence and the evacuated (i.e., coronal) part are taken homoge- 
neous, with densities pp and p^, respectively. The external, coronal medium has density pc, which 
is also homogeneous. Subscripts p, e, and c denote the prominence part, the evacuated region, and 
the corona, respectively. In general, the subscript denotes equilibrium quantities without refer- 
ring to a particular region. The density contrast of the prominence part with respect to the coronal 
plasma is a large parameter, withpp/pc = 200 a value usually considered. We assume that the evac- 
uated part of the tube has the same density as the corona. Hence, for a typical prominence density 
of Pp = 5 X 10"^^ kg m"^, the coronal and evacuated densities are Pc = Pe = 2.5 x 10"^^ kg m"^. 
In the prominence region, we include a transversely inhomogeneous transitional layer of thickness 
/, that continuously connects the internal prominence region to the external corona. A sinu soidal 
variation of the density is considered in the transitional layer (|Ruderman & Robe rtsI 120021). The 
limits l/a = and l/a = 2 correspond to a thread without transverse transitional layer and a fully 
inhomogeneous thread in the radial direction, respectively. The plasma in the prominence region 
is assumed partially ionized with an arbitrary ionization degree jUp. Both the evacuated part and the 
corona are taken fully ionized, hence /ie = fie =0.5. 

We use cylindrical coordinates, namely r, cp, and z for the radial, azimuthal, and longitudinal 
coordinates. The magnetic field is taken homogeneous and orientated along the z-direction, namely 
Bo = Boe^, with 5o = 5 G everywhere. The z-direction also coincides with t he axis of the cyli nder. 
For l/a = and full ionizatio n, the model is equivalent to that assumed by iDiaz et al.l (120021) and 
Dymova & RudermanI (l2005h . In these two works the prominence thread is located in the center 
of the cylinder. Here, we allow the thread to be displaced from the center of the tube. The length 
of the evacuated region on the left-hand side of the thread is L~, whereas the length of the right- 
hand side evacuated region is L^ . When L and Lp are fixed, we can express L^ = L - L~ - Lp, 
hence it is enough to select a value for L^ in order to set the length of both evacuated parts. 
The allowed values of L~ are in the range < L~ < L - Lp. For L~ = 0, the thread is totally 
displaced to the left-hand side end of the flux tube, while the contrary occurs for L" = L - Lp. 
For L~ = Lg = \[L- Lp) the prominen ce th read is located at the center of the tube, i.e., the 
configuration studied by iDiaz et al.l (|2002h and iDymova & RudermanI (|2005|) . For simplicity, we 
fix the origin of coordinates at the center of the prominence region, so that the interfaces between 
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the prominence plasma and the evacuated zones are located at z = +Lp/2. The photospheric walls 
are therefore located at z = z" „ = -Lp/2 - L~ and z = zl.u = Lp/2 + L^, with z. 



'■wall 



and z;,ii 

the position of the left and right walls, respectively. Thus, for Lp = L, i.e, a tube fully fill ed with 
prominence material, the model reduces to the configuration studied by ISoler et alJ (|2009q) . 



2.2. Basic Equations 



The governing MHD equations for a partially ionize d plasma are derive d in iForteza et al 



(120071) after the s i ngle-fluid tr eatment of iBraginskiil (| 19651) (see also details in lPintoet alj 12008 



Soler et all2009bl : ISolei]|2010() . Electrons, ions (i.e., protons), and neutral hydrogen are the species 
taken into account. We assume s mall-amplitude pe rturbations over the equilibrium state and 
the basic equations are linearized (ISoler et al.l l2009bl Equations (4)-(7)). In a partially ionized 
plasma, the induct ion equation cont ains diffusion terms related to the collisions between the dif- 
ferent species (see IPinto et al.ll2008b . For example. Ohm's diffusion is governed by electron-ion 
collisions, whereas Cowling's diffusion is do minated by ion-neu tral collisions. Here, we neglect 
Ohm's diffusion from the induction equation (ISoler et al.ll2009bl Equation (7)) because its role in 
a partially ionized prominence plasma is much less important than that of Cowling's diffusion. In 
addition, since we adopt the yS = approximation, gas pressure effects are neglected. 



Next, we follow iForteza et al.1 (|2008l) and lSoler et al.l (|2009bh and take a time dependence of 
the form exp (-icot), with co the oscillatory frequency, so that we can include the effect of Cowl- 



ing's diffusion in the definition of a modified Alfven speed squared as F^ = v^ - iioijco, where 



Ao 



■^Ao 



^Ao 



= B^/fxpo is the Alfven speed squared and t/q is the Cowling's diffusion coefficient, with 

fu = An X 10"^ N A~^. Th e expression of Co wlin g's diffusivity , rjc^, in terms of the equilibrium 

properties is given in, e.g.. ISoler et al.l (|2009bl) and lSoled (|2010|) . Thus, the relevant equations for 

our investigation are 

/Sv, 1 

(1) 



o 1 

Po^ = -(VxBi)xBo, 
ot n 



F^ 

4!^Vx(viXBo), 



^Bi 



V 



(2) 



Ao 



where Vi = (v,-, v^, vA and Bi = [Br, B^, B.) are the velocity and the magnetic field perturbations. 
Note that Vj. = in the^S = approximation, and Faq/vaq = 1 when Cowling's diffusion is absent. 

Equations ([T]) and ^ can be combined to arrive at the following equation for the total pressure 
perturbation, pj = BoB^/fi, namely 






rioV^PT = 0, 



(3) 



along with an equation relating the total pressure and radial velocity perturbations as 



df 



Ao 






1 5Vt 



(4) 



Po dtdr 

Note that Equation ([3]) is only valid in the regions with homogeneous densities, hence it cannot 
be applied in the transversely inhomogeneous transitional layer. Now, we write all perturbations 
proportional to exp {-icot + iniip), where m is an integer representing the azimuthal wavenumber. 
In the absence of magnetic twist, both positive and negative values of m are equivalent, so hereafter 
we restrict ourselves to positive values of m. For kink oscillations, m = 1. Equation ([3]) becomes 



d'^pj 



dpi 



1 d 

+ \r- 

dz r dr\ dr 



+ 



0) 

. Ao 



m 



2^ 



Pi = 0. 



(5) 



Since Pp, Pe, and p^ are uniform in their respective regions, the corresponding Cowling's 
diffusivities, ?7cp, rjcs, and rjcc, respectively, are also uniform. Both the corona and the evacuated 
region are fully ionized and much less dense than the prominence plasma, so we have rjcc «; ?7cp 
and T/ce -^ ?/cp- For the sake of simplicity, we set //ce = ?/cc = 0, and the effect of Cowling's 
diffusion is only considered in the prominence region. If considered. Cowling's diffusion in the 
evacuated and coronal regions would have a very minor influence since Cowling's diffusivities in 
these regions are much smaller than that in the prominence plasma. 



2.3. Mathematical Method 



For a fully ionized plasma, the general invest igation of t h e idea l transverse MHD oscillations 
supported by our equilibrium was performed by iDiaz et al.l (|2002h in the case II a = and for 
the prominence thread centered within the magnetic tube. These authors obtained the oscillatory 
frequencies and eigenfunct ions for arbitrary values of L, Lp, and a. Here, we could follow a 



treatment similar to that of iDiaz et al.l (|2002h . but this requires a significant mathematical effort 
beyond the pu rpose of the present investigat ion. Instead, we consider the much simpler approach 
introduced by iDymova & RudermanI (|2005|) . who studied the same configuration but in the thin 
tube (TT) limit, i.e., for alL «c 1 and a/Lp « 1. To check the validity of this approximation in 
the context of prominence thread oscillations, we take into account the values of a and Lp reported 
from the observations (see Sec. 12. II) and assume L ~ 10^ km. We obtain a/Lp and ajL in the ranges 
2 X 10"^ < a/Lp < 0.1 and 5 x 10 "^ < ajL < 3 x 10~^, meaning that th e T T approxiniation i s 
justified in pr ominence fine structure s . As s hown by iTerradas et al.l (|2008|) and Diaz et al.1 (120101) 



the me thod of iDymova & RudermanI (|2005l) shows an excellent agreement with that of 
(|2002|) when realistic values of a/Lp and alL are taken into account. 



Diaz et al. 
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Following iDymova & RudermanI (|2005|) . we can perfomi a different scaling of Equation ^ 
inside the tube and in the corona. For perturbations inside the tube, the characteristic scale in the 
r-direction is a, while the characteristic scale in the z-direction is L. Since a/L «c 1, the term with 
the longitudinal derivative and the term proportional to 0/ are much smaller than the other terms. 
In such a case, Equation ([5]) inside the tube reduces to 



dr\ dr I r^ 
with i = p or e. The solution of Equation ^ for regular perturbations at r = is 

PTi«Aia)|^-j , 

where Aj (z) is an arbitrary function of z- 



(6) 



(7) 



On the other hand, the characteristic scale of perturbations outside the tube, i.e., in the corona, 
is L in both the r- and z-directions, so that no terms can be neglected in Equation ([5]). However, we 
can express the total pressure perturbation in the corona as />tc = ^c iz) F (r) and use the technique 
of separation of variables to obtain the following expressions, 

d^F 



IdF 



dr^ r dr 



m 



+ z-n-\K + -T\F = o, 



and 



d^A, 









with 



Ae = 



at z = z: 



wall' 



(8) 



(9) 



where fc„ is a separation constant, with n an integer accounting for the different radial harmonics. 
In Equation @ we have taken into account that Cowling's diffusion is neglected in the corona. 
Equation ® is the modified Bessel Equation. Here, we only consider trapped 



sorL 



Ac- 



modes and assume kl > 0. This last condition may not be satisfied for high harmonics, i.e., 
large values of n, but we need not worry about this issue since here we focus our investigation 
on the fundamental mode, which is non-leaky in the present configuration. Then, the solution 
of Equation ([8]) is F (r) = K^ (Kr), with K^ the mo dified Bessel function of the second kind. 
An asymptotic expansion near the tube boundary (e.g. jAbramowitz & Stegunlll972l) allows us to 
express the total pressure perturbation in the corona as 



Ptc ~Ac(z)I-1 



(10) 



Next, our method follows closely that of lPymova & RudermanI (|2005L 120061) . For the sake of 
simplicity, we omit here the details, which are given in Appendix lAl After considering appropriate 
boundary conditions for the solutions of Equations ^ and (flOl) . the dispersion relation for kink 
oscillations damped by resonant absorption and Cowling's diffusion is obtained. 



3. DISPERSION RELATION AND APPROXIMATIONS 

The general dispersion relation for kink oscillations is 

g^^^osfe^)^"^fe-^e)-QeSin(^^)sin(^4-) 
^^^«in(a;^)'^o4^^e) + ^^eCOs(^^)sin(^L-) 

+ ^ ■ ^ ■ = 0, 

^^^sin(a^^)cosfe4-) + c.eCos(^^)sin(^4-) 



with 



Pp^'ap+P'^'ac _ Ij^^(r2 , „2 ^ ( PpPc: \ m/a 



'{[V^L)( 

\Pp+Pc / 



Si = -. ^^ (12) 



1 ' I Fppc 1 m/a 



-</ic(,^ 



t^Rl^rPolo 

m/a 



(11) 



^' = -^T-^ ^ —^ ^ ' (13) 

Pp^ Ap+P'^'ac • /p2 . ,.2 \ f PpPc \ m/a 
Pp+Pc "^ V Ap "^ ""Ac; I pp+p, I 6;R|fJ,po|„ 



2 Pp^ip+PcVL 2 PeVJe+PcVic ,... 

^ Pp+Pc Pe+Pc 

and h = ^ I - ^, with a>R the real part of the frequency and I^^PoL the radial derivative of the 

\ '^*p 

density profile evaluated at the resonance position, which has been approximated by the thread 
mean radius, a. The quantity ctp is here called the modified kink speed, which takes into account 
both the effect of Cowling's diffusion (through F^ ) and the effect of resonant absorption in the 
thin boundary (TB) approach. Extensive details are given in Appendix |Al If the terms related to 
resonant absorption are omitted, one has b^ = and c? becomes 

_2 Pprjp+PcvL ,,_, 

Cku = — ' (15) 



^^ Pp+Pc 



which reduces to the ideal kink speed, c? (Equation (fT4l)). when Cowling's diffusion is neglected, 

1-e-' rip = %• 

Equation (fTTl) is a transcendental equation that has to be solved numerically. Some analytical 
progress can be performed if the prominence thread is centered within the tube, i.e., L~ = L^ = 



10 



Ml - Lpj, and we focus on the fundamental kink mode. This solution corresponds to the mode 
with the lowest frequency. In such a case, Equation (fTT)) can be simplified to 



1 CjO 

tan 



cot 



C/te 



Cke 



L-U 



0. 



(16) 



The fundamental kink mode is given by the first root of Equation (1161) . A first-order Taylor ex- 
pansion for small arguments of the trigonometric functions of Equation (fT6l) provides us with an 
approximation to the frequency as 



( 



CO 



(l-l,)l, 



'-kp 






(17) 



We expect Equation (flTl) to be valid when both co and L^/L are small quantities, so that the argu- 
ments of the trigonometric functions of Equation (fT6l) remain small. Note that Equation (fTTl) fails 
to represent the kink mode frequency in the limits Lp/L ^ 1 and Lp/L -^ 0, so one should consider 
intermediate values of Lp/L in Equation (flTl) . The correct expressions for the fundamental kink 
mode frequency in these limits are 



OJ = —Ckp 



1 - -r, for Lp/L = 1, 



and 



TT 
(^ = jCke 



for 



Lp/L = 0. 



(18) 



(19) 



We can extract two main results from Equation (fT71) . First of all. Equation (fT71) only depends 
on the physical properties of the prominence region and the corona through c^, c^, and b, and 
includes no contributions from the evacuated part. And second, the form of Equation (fTT]) is similar 
to the approximation of the ideal kink mode frequency in a homogeneous tube, i.e., a/ « ^^c? , 
where k, is the longitudinal wavenumber. Thus, it seems that the main differences between the 
expression for the homogeneous tube and that for the partially filled tube are that 4/(L - Lp)Lp 
replaces k^, and that a redefined kink spee d has to be taken in t o acco unt. This approx i matio n of 
the frequency is similar to that obtained by Ijoarder & Roberts! (|1992h and lOliver et al.l (1 1993b for 



Diaz et al. 



( 20101) 



the string (or hybrid) modes of their slab configuration and to that obtained by 

in the context of thread seismology using period ratios. However, we must bear in mind that 

Equation (fTTl ) is only a first-order approximation to the kink mode frequency. 

Li the general case, i.e., when the prominence region is allowed to be at any position within 
the tube, one should consider the dispersion relation given by Equation (ITTI) . We can follow the 
same procedure as before and perform a first-order Taylor expansion for small arguments of the 
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trigonometric functions of Equation (fTTI) . Then, the following approximation for the frequency is 
obtained, 



4L 



(X> 



[(L-Lp)Lp + 4L-4-]Lp 



-^ 






(20) 



Note that the only information from the evacuated zones present in Equation (1201) is their lengths 
Lg and Lg , but no additional physical property of these zones contributes to Equation (|20l) . In the 
centered case, L~ = L^ = \^\L- L\ Equation (l20l) reverts to Equation (fTTI) . Now, we can consider 
the limits L^ ^ or L^T ^ 0, which correspond to the prominence thread totally displaced toward 
an end of the tube. In such limits. Equation (l20l) becomes 



4L 



iO 



(L-L,)Ll 



'/rp 






(21) 



The ratio of Equation (fTTI) to Equation (|2T)) estimates the shift of o/ when the thread is displaced 
from the center toward the end of the tube. Denoting this ratio as 6a/, we obtain 6a/ = Lp/L. 
Since Lp/L < 1, we expect the frequency of the kink mode to increase as the prominence thread is 
displaced form the central position. Obviously, for Lp/L = 1 there is no frequency shift because 
the prominence plasma occupies the whole tube. 



3.1. Period of the fundamental kink mode 



First, we focus our analytical investigation on the period. Here, we neglect the effect of the 
damping mechanisms since we assume that the kink mode period is only slightly affected by the 
presence of the damping mechanisms. As a first approximation, we consider Equation (l20l) with 
Ckp = Ckp and b = 0,so now a» is a real quantity. We compute the period as P = 2n/aj, obtaining 



p^Jl_ |Pp+Pc jK^-M^P+'^^e^el^P 



VAp V 2pp 



(22) 



Since the expression for the period is known, we could compare the theoretical periods with 
those observed and apply the MHD seismology technique to prominence fine structure oscilla- 
tions. However, in our case P depends on many parameters of the model, so Equation (|22)) alone 
is not very useful from a seismological point of view. If we assume that the prominence thread is 
located at the center of the magnetic tube. Equation (122)) becomes 



TT Pp + Pc 

VAp V 2pp 



V(^-^p)- 



(23) 
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In addition, in the case of prominences one has pp » Pc, so Equation (|23l) can be simplified to 



n 



V2vAp 



^(l-l,)l,. 



(24) 



Altho ugh the value of Lp can be measured from Ha observations of filaments (e.g. jLin et alj|2007l 
20091) . we still have two parameters, i.e., the total tube length, L, an d the prominence A lfven speed, 
that are both difficult to determine from the observations. Recently. iDiaz et al.l (|2010h showed that 
the ratio of periods of different overtones is a useful quantity to perform seismology of promi- 
nence threads, since additional parameters as, e.g., the prominence Alfven speed are dropped from 
the expressions (for details about the irn portance of the period ratio for coronal seismology, see 
the recent review by lAndries et al.ll2009h . Some additional remarks about the MHD seismology 
technique are given in Section [51 



3.2. Damping by Cowling's diffusion 



Here we study the kink mode damping. Let us consider the case without transverse transitional 
layer, i.e., l/a = 0, so the damping is exclusively due to Cowling's diffusion and the frequency is 
complex, a» = a»R + icoi, with ojr and oji the real and imaginary parts of the frequency, respectively. 
Then, c^ is given by Equation (fT5l) and b^ = since there is no resonant damping for l/a = 0. For 
simplicity, we consider that the prominence thread is located at the central position. Equation (fTTj ) 
allows us to obtain the real and imaginary parts of the frequency as 



0) 



- I 



Pp+Pc 
Pp'/Cp 



PpnCp 



1 



ll/2 



Pp+Pc/ (l-l^)l^ 



yl(L-L,)L, 



Pp+Pc/ (L-Lp)Lp 



(25) 



By setting the real part of Equation (|25l) equal to zero, we obtain two critical values of Lp/L, 
namely 



(^p/^)l 



1 

2 



1 - 



2pp 



Pp+Pc 



ll/2 



nk 



(26) 



with 77cp = Tjcp/vApL- Hence, the kink mode only exists for [Lp/Lj . < L^/L < [Lp/Lj . . We cast 
Equation (|26|) for pp/pc = 200 and the extreme case of an almost neutral plasma with jup = 0.99, 
obtaining f Lp/L j . « 10"^ and f Lp/L j . « 0.99999. For smaller values of /ip, (Lp/Lj . decreases 

and (Lp/L) . increases. Hence, the presence of these critical values is irrelevant for realistic values 
of Lp/L. 
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For Lp/L far from the critical values, one can drop the second term in the real part of Equa- 
tion (|25] ). We compute the damping time, td = l/|ct>i|, due to Cowling's diffusion as 



rD-2 



1 /Pp+Pi 



Pp'/Cp 



■)(^-^)^P 



(27) 



while the ratio of the damping time to the period is 



1/2 



Td _ J_ / Pp+Pc \ 1 
P 2n\ pp I fjcp 



ZE. 1 ZE 
L I L 



1211-^1^. 



(28) 



Forpp/pc = 200, Lp/L = 0.1, and L = 10^ km. Equation ([28]) gives td/P « 5 x 10^ for/ip = 0.8, 
and td/P « 150 for/ip = 0.99. These results indicate that, as in a homogeneous thread, an almost 
neutral prominence plasma is needed, i.e., /ip « 1, for the damping due to Cowling's diffusion to be 
efficient. Although the precise ionization degree is unknown, such large values of /ip are probably 
unrealistic in the context of prominences. 

It is straight-forward to extend Equation (1281) to the case in which the prominence region is 
not at the center of the tube, obtaining 



1/2 



Td _ 1 / Pp+Pc \ 1 
P 2n\ pp I ?7cp 



^p I ZE. + A^^^^ 
L L 



L2 



(29) 



Performing the limits L~ ^ or L^ ^ to this last expression, we can obtain the damping 
ratio when the thread is totally displaced toward the ends of the tube, namely (TD//')end- Then, 
we compare {Tu/P)end with the damping ratio obtained in the centered case from Equation (|28l) . 
namely (td//*) center- Thus, (td /P) center / (^d/^) end = -^Lp/L < 1. Therefore, the minimum value 
of the damping ratio by Cowling's diffusion takes place when prominence region is located at the 
center. 



3.3. Damping by resonant absorption 

Here, we study the general case l/a i^ 0. The full expressions of Ckp and b^ given by Equa- 
tions ([T2I) and (fT3]) are taken into account. Again, we assume that the prominence thread is centered 
within the magnetic tube. We use Equation (fTTl) to provide an expression for the ratio cui/ojr after 
neglecting terms ofOloj^j and Ofo^jL ^j. Thus, we obtain 



(Or 



7r(Pp-Pc) 



m/a 



Pv 



^ (pp+Pc) l^'-^ol" \Pp+Pc/ /2(i -b.\h 



1/2 



ricv 



(30) 



14 



The first term on the right-hand side of Equation (l30l) is caused by resonant absorption and the 
second term is due to Cowling's diffusion. The term related to Cowling's diff usion is also presen t 
in the case //a = 0. We can compare Equation (|30l ) with Equation (28) of ISoler et al.l (l2009c|) 
valid for a homogeneous tube. We see that both equations coincide if the replacement of k^ by 
2/ -yJ(L - Lp)Lp is done in their expression and t/cc = 0. Therefore, this fact suggests again that 
the results for the homogeneous tube can be extended to a partially filled tube by selecting the 
appropriate value for the longitudinal wavenumber. 

Note that the term related to the damping by resonant absorption in Equation (I3U1 ) takes the 
same form as in a homogeneous tube and does not depend on Lp/L. This result does not mean 
that the real and imaginary parts of the frequency do not depend on Lp/L, but both quantities are 
affected in the same wa y so that their ratio re main s unaffected. This irn portant result is consistent 
with the conclusions of lAndries et al.l (|2005h and lArregui et al.l (|2005h . who found that the kink 
mode frequency of a longitudinally stratified tube is the same as that obtained for a homogeneous 
tube with density Pmean, with Pmean the mean density weighted with the wave energy. Since in 
our equilibrium the transversely transitional layer is only present in the dense part of the tube and 
both the evacuated zone and the corona have the same density, resonant absorption only takes 
place in the dense part of the tube. Therefore, the mean density of the pa rt of the tube where 
reso nant absorption take s place is, obviously, p^ean = Pp- Hence, according to lAndries et al.l(|2005|) 
and lArregui et al.l (|2005|) . the kink mode damping ratio in our case must be the same as that of a 
homogeneo us tube with density Pp, as our results indicate. This conclusion is also equivalent to that 
obtained by lOymova & RudermanI (|2006r) . who showed that the damping ratio in a longitudinally 
inhomogeneous tube in the TT and TB approximations does not depend on the particular form of 
the longitudinal density profile if the density contrast between the internal and external plasmas is 
constant. 



By assuming a sinusoidal density variation in the transitional layer, the expression for td/P 
according to Equation (I3U1 ) is 



P 



1-1 



in—- 



Pp 



aPp+Pc \Pp+Pcl /2(i -b.]h 



1/2 



(31) 



which is equivalent to Equation (29) of ISoler et al.l (|2009d) for fjcc = and k^ = 2/ ^{L - Lp)Lp. To 
perform a simple application, we compute r^/P from Equation (|3TI) in the case m = 1, Lp/L = 0.1, 
L = 10^ m, and l/a = 0.2, resulting in td/P ~ 3.18 for a fully ionized thread (flp = 0.5), and 
td/P ~ 3.16 for an almost neutral thread (fip = 0.95). We n ote that the obtaine d damping times are 
consistent with the observations. Moreover, as obtained by ISoler et al.l (|2009ch . the contribution of 
resonant absorption to the damping is much more important than that of Cowling's diffusion, so 
the ratio td/P depends only very slightly on the ionization degree. Then, the second term on the 



-15- 

right-hand side of Equation (|3TI) can be neglected and Equation (|3TI) becomes 

td _ 2 1 1 Pp + Pc 
P nml/app - Pc' 



(32) 



which coincides with the expressions provided by lRuderman & RobertsI (|2002h and lGoossens et al. 



(120021) in the case of longitudinally homogeneous coronal loops. Equation (l32l) can be further 



simplified by restricting ourselves to the kink mode (m = 1) and for pp :» pc, hence 

P nil a 

meaning that the transverse inhomogeneity spatial- scale can be estimated if both the period and 
damping time are measured from the observations. By combining Equations (l24l) and (1331) . we get 
the expression of the damping time by resonant absorption as 

r„.^2;V(iri:)i;. (34) 

Finally, we can compute the damping ratio when the prominence thread is displaced from the 
center of the tube. In such a case, only the term related to Cowling's diffusion is modified in the 
way shown in Section [3^ while the term related to the resonant damping is not affected at all. The 
damping time by resonant absorption in the general case is 



V2 1 [(L-Lp)Lp + 4L-L,^]Lp 
VAp II a M L 

which has the same dependence on L^ and L^ as the period (Equation (1221)'). Since the resonant 
damping dominates over Cowling's diffusion, we anticipate by means of this analytical estimations 
that the damping ratio is almost unaffected by the position of the prominence region within the fine 
structure. 



4. NUMERICAL RESULTS 

4.1. Centered prominence thread 

In this section, we assume that the prominence region is centered within the tube, i.e., L^ = 
Lg = \\L- Lp j. We numerically solve the dispersion relation (Equation (fT6l) ) by means of stan- 
dard methods and obtain the frequency of the fundamental mode. We study the dependence of the 
results on Lp/L. 
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0.2 0.4 



Fig. 2. — A{z) function (in arbitrary units) corresponding to the fundamental kink mode, (a) Results 
for Lp/L = 0.1 (solid), 0.5 (dashed), and 0.9 (dotted) when the prominence thread is located at the 
central part of the magnetic tube, (b) Results with Lp/L = 0.2 for L~IL = 0.2 (dashed), 0.4 (solid), 
and 0.7 (dotted). The thick part of the lines in panel (b) denotes the position of the prominence 
thread. 

We plot in Figure [2^ the A{z) function corresponding to the fundamental even kink mode for 
different values of Lp/L. The A{z) function gives the dependence of the transverse displacement 
in the longitudinal direction at r = a. We see that A{z) is mainly confined within the prominence 
part of the flux tube and satisfies the line-tying condition at z = z*^,j = +L/2. For a homogeneous 
prominence tube, i.e., Lp/L = 1, the A{z) function becomes proportional to cos f|zj. 



4.1.1. Case without transverse transitional layer (I I a = 0) 



First, we take into account the case without transverse transitional layer, i.e., l/a = 0, and so 
we study the kink mode damping due to Cowling's diffusion exclusively. Figure [3^ displays the 
period, P, as a function of Lp/L for different values of the ionization degree in the prominence 
region, whereas Figure [3b shows the corresponding values of the damping time, td. Both values 
are given in dimensionless form with respect to the internal Alfven travel time, tap = L/vAp. We 
see that P/tap increases as Lp/L increases, and tends to the value for a homogeneous prominence 
cylinder when Lp/L ^ 1 . In addition, we find that P is independent of the ionization degree. On the 
contrary, td is strongly dependent on the ionization degree, as expected, td/tap slightly increases 
as Lp/L becomes larger. We see that the analytical expressions for the period (Equation (l24l) ) and 
the damping time (Equation (ITTI) ) are in agreement with the solution of the full dispersion relation 
for realistic, small values of Lp/L, i.e., Lp/L < 0.4, whereas the approximate expressions diverge 
from the actual solution when the prominence region occupies most of the magnetic tube. As 
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Fig. 3. — Results in the case without transverse transitional layer and for the prominence thread 
located at the central part of the magnetic tube, (a) Period, P, of the fundamental kink mode in 
units of the internal Alfven travel time, TAp, as a function of Lp/L. The horizontal dotted line 
corresponds to the period of the kink mode in a homogeneous prominence cylinder. The symbols 
are the approximation given by Equation ((24)) . (b) Damping time, td, in units of the internal Alfven 
travel time, TAp, as a function of Lp/L. The different lines denote /ip = 0.5 (dotted), 0.6 (dashed), 
0.8 (solid) and 0.95 (dash-dotted). The symbols are the approximation given by Equation (ITT]) for 
/ip = 0.8. (c) Td//' versus Lp/L. The line styles have the same meaning as in panel (b), and the 
symbols are the approximation given by Equation (1281) . 
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Fig. 4. — Results in the case with a transverse transitional layer and for the prominence thread 
located at the central part of the magnetic tube, (a) td, in units of the internal Alfven travel time, 
TAp, and (b) td/P as a functions of Lp/L. The different lines in both panels denote Ija = 0.05 
(dotted), 0.1 (dashed), 0.2 (solid) and 0.4 (dash-dotted). The symbols in panels (a) and (b) are the 
approximations given by Equations (l34l) and (l32l) . respectively, with Ija = 0.2. 



commented before, high-resolution observations suggest that the parameter Lp/L is small in the 
fine structures of prominences. On the other hand. Figure [3]: displays Td/L* versus Lp/L. The 
numerical solution of the dispersion relation shows little dependence on Lp/L, while the analytical 
approximation (Equation (|28l) ) diverges from the numerical value in the limit of large Lp/L. Given 
the large values of td/P obtained, we can conclude that the efficiency of the damping due to 
Cowling's diff"usion in a pa rtially filled flux tub e does not improve with respect to the longitudinally 
homogeneous tube case of ISoler et al.l (|2009blJcl) . 



4.1.2. Case with a transverse transitional layer (Ija i^ Oj 



Now, we take the case Ija 4^ into account. The kink mode is damped by resonant absorption 
in the transverse transitional layer. We have computed both the period and the damping time of 
the fundamental kink mode as a function of the diff'erent parameters, namely /ip, //a, and Z.p/L. 
Regarding the period, we find that both its value and its dependence on Lp/L are the same plotted 
in Figure[3ti because the period is almost independent of /ip and Ija. For the sake of simplicity, we 
do not repeat this Figure again and refer to Figure [3^. Therefore, the presence of the transverse 
transitional layer in the TB approximation does not modify the period of kink oscillations with 
respect to the case II a = 0. The period could be slightly affected if thick layers, i.e, l/a > 1, are 
considered and the resistive equations are solved numerically instead of assuming the TB approach 
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see, e .g., I Van Doorsselaere et al.l |2004|) . This issue will be further addressed by lArregui et al 



2OIOI) . 



Figure H^ shows Td/tap versus L^/L for different values of l/a. These computations corre- 
spond to an ionization degree /ip = 0.8, but equivalent computations for other values of /ip provide 
almost identical results because the effect of Cowling's diffusion is negligible in comparison to 
that of resonant absorption. As expected, the value of the damping time decreases with //a. The 
approximate value of tq given by Equation (|34l ) is in good agreement with the full solution for 
Lp/L < 0.4, as happens for the period. In order to assess the efficiency of the resonant damp- 
ing. Figure lib displays the corresponding values of t-qIP. In comparison to the damping ratio by 
Cowling's diffusion (see Fig. [3t), much smaller values of t-qIP are now obtained. As predicted 
analytically by Equation (|32l ), Td//' is almost independent of Lp/L. By comparing Figures |3^ and 
S^, we see that both the period and the damping time have a very similar dependence on Lp/L, so 
the dependence on Lp/L is canceled when the damping ratio is computed. In this case, a very good 
agreement between the numerical result and the analytical approximation (Equation (1321) ) is found 
even for large values of Lp/L. 

4.2. Effect of the position of the prominence thread within the magnetic tube 

In this Section we study the effect of the position of the prominence region within the magnetic 
flux tube. The results of the previous Section 14.11 correspond to the case in which the thread 
is located at the center of the cylinder. Here, we allow the dense region to be displaced from the 
center of the tube. Hence, we must consider the general dispersion relation given by Equation (fTTI) . 
The dispersion relation is solved numerically for the lowest frequency solution, equivalent to the 
fundamental kink mode of the centered case. 

We display in Figure |2j5 the A{z) function for different values of L~IL for Lp/L = 0.2 and 
/ip = 0.8. Since the oscillation is dominated by the prominence physical properties, we see that 
the maximum of A(z) is always in the prominence region (denoted by the thick part of the lines in 
Fig.[2j5), regardless of its location within the flux tube. 

Next, we plot the period (Fig. [5^) and the damping time (Fig. [Sh) as functions of L^jL for 
Lp/L = 0.1, /ip = 0.8, and different values of 1 1 a. We obtain that the longest period takes place 
whentheprominencethreadiscentered within the flux tube, i.e., when Lg/L = ^(l - Lp/Lj = 0.45 
for this particular set of parameters, and P decreases symmetrically around L^/L = 0.45 when 
L~/L increases or decreases. The dependence of r^ on L~ /L shows the same behavior as P. 
Such as happens with the dependence on Lp/L (Fig. H]), the dependence on L~ /L also cancels 
out when the damping ratio is computed (see Fig. [5t). Hence, in our model the value of Tu/P is 
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Fig. 5. — (a) Period, P, of the fundamental kink mode in units of the internal Alfven travel time, 
TAp, as a function of L'jL. The symbols are the approximation given by Equation (|22l) . (b) 
Damping time, tq, in units of the internal Alfven travel time, TAp, as a function of L^/L. The 
different lines denote Ija = 0.05 (dotted), 0.1 (dashed), 0.2 (solid) and 0.4 (dash-dotted). The 
symbols are the approximation given by Equation (l35l) with l/a = 0.2. (c) td/P versus L^/L. The 
line styles have the same meaning as in panel (b), and the symbols are the approximation given by 
Equation (|32l) with l/a = 0.2. In all computations, Lp/L = 0.1. 
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independent of both Lp/L and L'/L. We can also see in Figure [5]: that the approximate td/P given 
by Equation (|32l ) remains valid even when the thread is not located at the center of the magnetic 
tube. 



5. IMPLICATIONS FOR PROMINENCE SEISMOLOGY 



The obtained results have direct implications for the determination of physical parameters in 
prominence fine structures using MHD seismology together with observed periods and damping 
times. The technique of MHD seismology has been previously applied by some authors to ob- 
tain information of the plasma physical conditions in the context of coronal lo op oscillations (e.g., 
Nakariakov & Ofmaj200ll:lArregui et allboOTl boOSal: lOoossens et allboosl) . prominence global 



Robertslll991 



oscillations (e.g., 

oscillations (e.g. jTerradas et al 



Regnier et al.l 



2008 



Lin et al. 



200ll:|Pouget et al.l 120061) . and prominence thread 



20091). As partial ionization has a negligible effect 



on the damping of oscillations we here concentrate on the inversion of parameters using theoret- 
ical results for resonantly damped ei genmodes. Following the analytical and numerical inversion 
schemes by lArregui et al.l (120071) and lGoossens et al.l (l2008b for coronal loops. lArregui & Ballestei 
(120 10|) have recently presented the inversion of prominence Alfven speed, transverse inhomogene- 
ity length-scale, and density contrast in oscillating filament threads using results for resonantly 
damped kink oscillations in one-dimensional (ID) filament thread models. A similar procedure 
can be followed with two-dimensional (2D) threads by considering the impact that the length of 
the thread has on the period and damping ratio of standing kink modes in partially filled fine struc- 
ture oscillations. 

In the following computations, we consider that the prominence thread is located at the center 
of the magnetic flux tube. In the 2D case and for small Lp/L, i.e, Lp/L < 0.4, the oscillatory period 
is approximately given by Equation (1231 ) for arbitrary density contrast, and by Equation (I2?l) in the 
limit of large density contrast. In the ID case, the kink mode period in the TT approximation is 



(36) 



where A is the wavelength. By comparing Equations (|23]) and (l36l) . we see that the expression in 
the ID case is equivalent to that in the 2D case if a particular or effective value of A is considered, 

namely A x n J(L - LplLp. This effective value of A allows us to generalize the ID inversion to 
a 2D configuration. Equation (l23l) is not accurate enough when Lp/L > 0.4, but we can still use a 
similar expression for the period for any Lp/L as 
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Fig. 6. — (a) Inversion of physical parameters in the ip^lpc, 1 1 a, VAp) space for a prominence thread 
oscillation with P = 20 min and td/P = 3, and for different values of Lp/L (indicated within 
the Figure). The thin continuous, dotted, and dashed lines correspond to the projections of the 
three-dimensional curves to the (Pp/pc, l/a)-, (l/a, VAp)-, and (vAp, Pp/pc)-planes, respectively, (b) 
Inversion of the prominence Alfven speed, VAp, as a function of the period, P, in the limit of large 
density contrast. The different lines correspond to Lp/L = 0.2 (solid), 0.4 (dotted), 0.6 (dashed), 
and 1 (dot-dashed). The shaded zone corresponds to the range of typically observed periods in 
thread oscillations, (c) Magnetic field strength, Bq, as a function of the prominence thread density, 
Pp, assuming an oscillatory period of P = 26 min. The different lines have the same meaning as in 
panel (b). In all these computations, L = 10^ km. 
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with C a correction factor that is computed by comparing the period given by Equation (1231) to 
the general result of Figure [3^ (solid line), obtained by solving the full dispersion relation (Equa- 
tion (fT6l)). Thus, A X Cn J(l - Lpj Lp in the general case. For example, by assuming L = 10^ km, 
we obtain C « 1.05 and A ~ 9.89x 10^ kmfor Lp/L = 0.1, whereas C ~ 1.58 and A « 1.98x 10^ km 



for Lp/L = 0.8. As expected, A ^> 2L = 2x 10 km in the limit Lp/L -^ 1. 

On the other hand, the damping ratio is unaffected by the length of the thread and the same 
expression for ID models given by Equation (l32l) holds. By following the analytical inversion 
scheme by iGoossens et al.l (|2008|) the valid equilibrium models that explain equally well a given 
set of parameters (P, Td, and L) are obtained. The resulting one-dimensional curve in the three- 
dimensional parameter space is shown in Figure [6^, for different values of the length of the thread. 
Regardless of the value of Lp/L, the inversion curve allows us to obtain well constrained values for 
VAp and l/a in the limit of high density contrast values. Because of the decrease in A produced by 
the decrease of the length of the thread, the obtained Alfven speed in the prominence decreases as 
Lp/L gets smaller. An accurate estimate of the length of the thread, in comparison to the length of 
the magnetic flux tube, is therefore crucial for the determination of the Alfven speed in the thread. 
As pointed out by iDiaz et al.l (|2010|) . Lp/L can be obtained from the ratio of the fundamental mode 
period to that of the first harmonic, if these values are reported from the observations. However, 
because of the independence of the damping ratio on Lp/L, the projection of the solution curve 
onto the (pp/pc, //a)-plane remains unaltered regardless of the value of Lp/L. Hence, the inverted 
value of the transverse inhomogeneity length-scale, which is obtained taking the limit pp/pc -^ oo, 
is not affected by different values of the length of the thread. 

The value of P has a direct impact on the seismological determination of the Alfven speed, 
with the inversion of Figure [6^ corresponding to P = 20 min. The effect of the period on the 
determination of the prominence Alfven speed is indicated in Figure |6b, where we see that the 
Alfven speed decreases as the period gro ws. Here, it is impor t ant recalling that the usually re 



porte d periods of oscillating threads (e.g., iLin et al.ll2007L 120091: lOkamoto et al.ll2007l : iNing et al 



20091) are in the range 2-10 min. Note the very large values of VAp obtained in the range of 
observed periods (shaded zone in Fig. |6b). By assuming pp = 5 x 10"'^ kg m~^ and 5o = 5 G 
as typical values of the density and magnetic field strength of quiescent prominences that can be 
found in the literature, respectively, the corresponding Alfven speed is VAp ~ 63 km/s, which for 
L = 10^ km and Lp/L = 0.2 gives P « 23.5 min. Hence, to obtain realistic values of VAp, we have to 
consider larger periods than those usually observed. Equation (1221) indicates that, for the same set 
of parameters, the period decreases if the thread is displaced from the center of the magnetic tube 
(see Fig.|5^). The maximum shift of the period with respect to the value in the centered case is pro- 
portional to -yjLp/L. Adopting the same parameters as before, the period varies from P « 23.5 min 
to P « 9.4 min when the thread is displaced from the center to the end of the magnetic tube, and 
so the period enters within the observed range. However, there is no solid basis to assume that all 
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short-period oscillating threads are located at the ends of their magnetic tubes. Alternatively, the 
presence of flows can also shift the oscillatory period. Although m ass flow has no t been included 
in our model, the time-dependent simulations of flowing threads by lTerradas et al.l (|2008|) indicate 
that flow has a minor influence on the period because the flow velocities are much smaller than the 
Alfven speed. In active region prominences, the larger magnetic field strength could cause larger 
Alfven speeds, hence thread standing oscillations in active region prominences may have shorter 
periods than in quiescent prominences. 



Moreover. iLin et all (|2007ll2009h and lNing et all (|2009|) reported that the wavelengths of these 
short-period oscillations are in the range 700 - 8,000 km. The observed wavelengths are between 
1 and 2 orders of magnitude smaller than the wavelengths corresponding to standing oscillations 
computed from Equations (|36l ) and (|37l ). Thus, the reported short wavelengths are impossible to 
reconcile with the fundamental standing mode of the fine structure. These results suggest that the 
observed short-period and short-wavelength oscillations of threads in quiescent prominences may 
not be consistent with an interpretation in terms of standing kink oscillations. A more likely expla- 
nation of these short- p eriod and s hort-w avelength oscillations in terms of propagating kink waves 
has been performed by lLin et al.l (|2009|) . These authors seismologically inferred realistic values of 
the Alfve n speed and ma gnetic field strength by assuming a propagating wave interpretation (see 
details in lun et al.ll2009 ). 



The limited duration of the currently available Doppler time series may prevent the observa- 
tion of standing m odes with period s larger than 10 min (e.g., the time series last for only 18 min in 
the observations of lLin et al.ll2007|) . However, there are a few evidences of larger periods in longer 
time series that could be consistent with standing modes. These longer periods have been obtained 
from Doppler signals that have been averaged over a large area. As the spatial scales of standing 
modes are very large, their oscillatory patterns could be found in spatially averaged signals. On the 
contrary, the averaging process could mix signals coming from difl^erent adjacent thread s, and so 



there is no full confidence that the period corresponds to an individual thread oscillation. lYi et al. 



(| 19911) detected thread oscillations with period of 16 min and minimu m wavelen gth of 2 x 10"^ km 
in their Doppler observations with low spatial resolution (1"), while Linl (|2004() reported 26 min 
and minimum wavelength of 4 x 10"* km in their averaged Doppler signals. Although not only the 
periods but also the wavelengths reported in these two works are consistent with a standing oscil- 
lation, we have to be very cautious with this interpretation. Additional information as, e.g., the 
polarisation of the oscillations and the phase diff"erence between perturbations, would be needed 
for a more robust analysis and the unequivocal determination of the wave mode. 



To perform a simple application, let us assume that the period of 26 min reported by Lin 



(|2004|) corresponds to a standing thread oscillation. Then, the estimation of th e magnetic field 
strength in the prominence thread is possible. Following the analysis presented by lLin et al.l(|2009|) 



25 



using the observational period and the seismologically determined Alfven speed (as in Figl6h), the 
magnetic field strength can be computed for a given prominence density. This result is plotted in 
Figure [6]:. For a typical density of Pp = 5 x 10"'^ kg m""*, magnetic field strengths in the range 
4 - 7 G are obtained, approximately, when varying Lp/L between 0.2 and 1. These values of the 
magnetic field strength are in agreement with previous mag netic field measurements in quiescent 
prominences using the Hanle effect (e.g. JLeroy et al.lll984|) . Hence, the method can be applied in 
the future using real data if reliable observations of standing thread oscillations are reported. 



6. CONCLUSION 



In this paper, we have investigated standing kink oscillations of prominence fine structures. 
The longitudinal nonuniformity has been taken into account by modeling the fine structure as a 
magnetic tube only partially filled with the prominence material. We have followed an analytical 
method based on the TT approximation and have found a dispersion relation for kink oscillations 
damped by Cowling's diffusion and resonant absorption in the TB approach. This dispersion re- 
lation has been numerically solved and a parametric study of the solution has been performed. In 
addition to the general dispersion relation, we have obtained simple analytical approximations to 
the period, the damping time, and their ratio. 

Both approximate and full results conclude that resonant absorption is much more efficient 
than Cowling's diffusion for the kink mode damping, with the values of td/P in agreement with 
those reported in the obs ervations. As happe ns for long- wavelength propagating waves in the dense 
part of the fine structure (ISoler et al.ll2009a) . the prominence plasma ionization degree turns out to 
be irrelevant for the resonant damping of the oscillations. In addition, the value of T^tlP is found 
to be independent of both the position of the prominence thread within the magnetic tube and the 
length of the pron iinence region, and coincides with the va lue for a homogeneous and fully ionized 
prominence tube (lArregui et al.ll2008al : ISoler et al.ll2009a|) . 



Finally, we have discussed the seismological implications of our analytical results, in par- 
ticular Equations (l24l) and (|33] ). With these expressions, it is possible to estimate some relevant 
physical parameters of oscillating threads if the values of the period and damping time are avail- 
able from the observations. Following this idea, we have performed a seismological inversion of 
the prominence thread Alfven speed and the transverse inhomogeneity length-scale by using our 
theoretical results and adopting ad-hoc values for the period and damping time. We have shown 
that for short-period (2 - 10 min) and short-wavelength (700 - 8,000 km) thread oscillations, the 
determined Alfven speeds are much larger than the expected, realistic values, pointing out that 
short-period and short-wavelength thread oscillations may not be consistent with a standing kink 
mode interpretation and could be related to propagating waves. On the contrary, thread oscillations 
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with periods larger than 10 min and wavelengths larger than 10"^ km may be interpreted as stand- 
ing oscillations. In this last case, the Alfven speed and magnetic field strength estimated by the 
seismological inversion are realistic in the context of prominences. Thus, the method can be put 
into practice to extract indirect information of prominences when standing thread oscillations are 
unequivocally observed and the oscillation parameters, i.e., period, damping time, and wavelength, 
along with the thread length are provided from the observations. 

In this work, we have assumed that there is an abrupt jump of the density at the boundary 
between the prominence thread and the evacuated part of the magnetic tube. This simplifica- 
tion has allowed us to proceed analytically. Actually, one should expect a continuous variation 
of the plasma properties in the longitudinal direction between both regions, which could afl'ect 
somehow our present results. The study of the damping of kink oscillations in fu lly nonuniform 



two-dimensional fine structures is broached numerically by lArregui et al.l (120 10|) in a following 
investigation. Other additional ingredients as, e.g., magnetic twist or curvature and the presence of 
flows might be included in future works. 
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A. DERIVATION OF THE DISPERSION RELATION 

Here, we give extensive details about the method that leads us to the dispersion relation (Equa- 
tion (fTTT)). 



A.l. Boundary conditions atr = a 

First, we must consider appropriate boundary conditions for the solutions of Equations © 
and (fTUl ) at the cylinder edge, i.e., r = a. In the evacuated part of the tube, we assume Pe = Pc 
and there is no transverse tra nsitional layer. Hence, the boundary conditions are those given by 



Dymova & RudermanI (|2005h in their Equation (4), namely 



[bT]]=0, [[v,]]=0, at r = a for |z| > Lp/2, (Al) 
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where [[X]] stands for the jump of the quantity X. 

On the other hand, in the prominence part of the tube we consider the effec t of resonant ab- 
sorption in the transitional layer. We follow the treatment bv lAndries et al.l (l2005h. wh o generalize 
the concept of the jump conditions at t he resonance position of lSakuraiet alj (| 199 lb to the case 
of a longitudinally inhomogeneous tube. lAndries et al.l (|2005h combined the jump conditions with 
the TB approximation, i.e, l/a <^ 1, to obtain analytical expressions of the dispersion relation and 
the frequency for lo ngitudinally stratified tubes. The accuracy of this analytical method was nu- 
m erically verifie d by lArregui et al.l (l2005b . who found a good agreement between the expressions 
of lAndries et al.l (l2005h and the i r nurn erical computations. In the thin tube approximation, i.e., 
a/L «c 1 . iDymova & RudermanI (|2006|) follow a similar formalism and also provide equivalent ex- 
pressions for the jump conditions that can be applied to our perturbations. Hence, it is convenient 
to express pj and v^ as 

oo oo 

where pjn and v™ are the coefficients of the series expansions of pj and v,-, respectively, with 
respect to the functions G„ determined by the Sturm-Liouville problem 

d-G„ 



vi«,,. 



-AI (r) G„, 



(A3) 



with appropriate boundary conditions for G„ at z = ±Lp/2, with Al the corresponding eigenvalues. 
Equation (|A3I) describe the spectrum of Alfven modes, with An (r) the corresponding frequencies 
of the Alfven continuum. In general, it is not straight-forward to deduce the boundary conditions 
for G„ at z = +Lp/2 because they are given by the continuity of G„ at z = +Lp/2, and the value of 
G„ at z = +Lp/2 is also determined by the properties of the evacuated region. In a longitudinally 
homogeneous tube, i.e., for Lp = L and with the Alfven speed depending on the radial direction 
only, we simply have that the boundary conditions are G„ (+L/2) = and obtain An (r) = coa (r) = 
^VA(r),withn = l,2,... 



In our notation, the jump conditions for pjn and v™ provided by lDymova & RudermanI (|2006|) 



are 



tn /a 
[[PTn]] = 0, [[Vrn]] = -ncoRj—-;-^ — pxn, at r = r^n for \z\ < Lpll, 



(A4) 



■|PoA„l,^„ 

where r^n is the Alfven resonance position for the nth mode and A„ = ^ fa)| - /l^j, with ojr the 
real part of the frequency. According to Equation (|A2|) . the condition for p^n in Equation (|A4|) 
leads to [[px]] = at r = r^n for \z\ < Lp/2. On the contrary, the condition for v„ depends on 
|poA„|,.^, and a more detailed analysis is needed. 

For our subsequent analysis, we do not need the precise value of An (r) but only its functional 
dependence on the radial direction. For the given sinusoidal density profile in the transitional layer. 



28 



we can express the Alfven speed squared in the transitional layer as vi (r) = v\ If (r), with 



"Ap/ 



^«4{('%^ 



1 - ^ I sin 



j(r-a) 



(A5) 



Hence, Equation (IA3h is rewritten as 



V 



Ap- 



d_^ 
dz2 



= -i^(r)/(r)G„. 



(A6) 



With no loss of generality, we can assume that G„ is only a function of z, i.e., the different magnetic 
surfaces are not coupled to each other. So, according to Equation (IA6I ), the quantity Al (r) f (r) 
corresponds to the Alfven eigenvalue squared in the prominence part of the tube. Since Pp is 
homogeneous, its corresponding Alfven eigenvalue does not depend on r, meaning that the radial 
contribution of Al (r) and /(r) cancel out. Thus, we define A^^ = Al (r)f(r), with Ap„ a constant 
corresponding to the Alfven eigenvalue in the homogeneous prominence thread. Therefore, we 
have 

^^(r) = ^, (A7) 

where all the radial dependence of Al (r) comes from the function /(r). With the help of Equa- 
tion (IA7I) . we obtain that A„ = A^ (r) f (r) //(r), where the prime denotes the radial derivative. 
Finally, we use the resonant condition, namely A^ {r^n) = <^r, and write 



IPoA„ 



.A„ = ^R \drP^\r^, 



(A8) 



Thus, the condition for Vm in Equation (IA4I) becomes 



[[v™]] = -n — ■ — pTn, at r = r^n for \z\ < Lp/2. 



(A9) 



The value of rA„ is in principle different for each value of n and could be determined from the 
resonant condition A?^ {rp^n) = <^r if the eigenvalues A^ (r) were a priori known and the number of 
Alfven eigenmodes that are resonant to the kink mode is also known. Hence, the derivative of the 
density profile at each resonant position could be computed. A reasonable assumption in the TT 
and TB limits is to consider ta,, « a for all n, so that \d,-pQ\^ is a constant independent of n, meaning 
that we are assuming that all resonances take place at the same position. For our sinusoidal profile, 
\drPo\a * ^(Pp ~ Pc] /2/. Note that this is not a very strong restriction for the fundamental kink 
mode, since it is likely that the resonant condition is satisfied for n = I only, because A^ (r) grows 
as n increases and so only one resonance takes place. The approximation ta,, « a might not be valid 
for the kink mode overtones, but here we restict o urselves to the fundam ental mode. Alternatively, 
a simpler linear density profile could be adopted (jGoossens et al.ll2002r) in which the derivative of 
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the density profile does not depend on the position. Therefore and using Equation (IA2I) . we arrive 
at 

2 / 2 

m la 

' OjR\drPQ\J 



WvrW = -n , ,;, ^ , Pt, at r = a for |z| < Lp/2. (AlO) 



ForLp = L the jump condition of Equation (lAlOD consistingly reduces to that provided by lSakurai et al 



( 1991 ) 



A.2. Solution in the evacuated regions 

Let us consider first the boun dary conditions for the evacu ated parts (Equation (|A1I) ). The 



analysis here is identical to that of iDymova & RudermanI (|2005h . For the condition on the total 



pressure perturbation, we obtain Ag (z) = Ac (z) = A (z). Next, we rewrite Equation dH) as 

oz^ Pq or 

We evaluate Equation (lAllI) forr^a on both sides of the tube boundary. Thus, in the evacuated 
part. 



<xf + ^^^ = — ^- (All) 



^Ae^^ + ^ I're = ^ (z) , for r ^ Q, (A12) 



whereas in the corona. 



d^Vrr T iCt» m 
vie— f + a;Vc = A(z), for r ^ a. (A13) 

OZ^ Pc « 

According to the boundary condition for v^ given by Equation (|A1I) . Vre = Vrc. Thus, we combine 
Equations (IA12I) and (IA13I) to find the following two expressions 



5V . m p pe+pc 

til \ 2 

v. = -/ A(z), (A15) 

a cope -pc 

where we have considered that the magnetic field is homogeneous. Now, we diff'erentiate Equa- 
tion (IA15I ) with respect to z twice and compare the resulting expression with Equation (IA14I) . We 
obtain 

^ + ^Ab)=0. (A16) 

The quantity c^ corresponds to the kink speed in the evacuated region (Equation (fT4l)). Note that 
in our particular application pe = Pc, so c^ = Vae- However, we keep the general notation Cke in the 
following expressions. 
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To solve Equation (IA16I) . we consider the line-tying condition at the photosphere, i.e., A (z^^^A 
0. Therefore, the solution in the two evacuated zones is 



A(z) = < 
where Ci and C2 are constants 



' Ci sin [^ (z - Lp/2 - L:)] , for z > Lp/2, 
_ C2 sin [^ (z + Lp/2 + 4-)] , for z < -Lp/2, 



(A17) 



A.3. Solution in the prominence thread 

In the prominence thread, we adopt the TB approach and use the jump conditions given by 
Equation (IA4I) as our boundary conditions. The combinatio n of both formalisms to investigate 
resonant waves in coronal flux tubes has been reviewed by iGoossens et al.l (|2006|) . Again, the 
condition over the total pressure perturbation gives Ap (z) = Ac (z) = A (z). Near the boundary 
we express Vyc = v^p + 6vr, where 6vr is the jump of the radial velocity perturbation provided by 

Equation (|A10I) . namely 

2 / 2 
ma 

-Pi. (A18) 



6Vr 



-n- 



As before, we evaluate Equation (|4]) on both sides of the tube boundary and, after some algebra, 
we arrive at the following expressions 



^ rp 



m 



10) Pp + Pc 



dz^ 



+ 



a V 

m^ la 



Ac 
2 /^2 



rip PpPc 

n 



A{z) 



OJR \drP^\a Vl, - E^p 



'Ac 



d'A (z) 
dz2 



+ co^A (z) 



(A19) 



'rp 



m i Pp^Ap + Pci^Ac 

"^^PpPc(vL-rip) 



ni^/a^ 



A(z) 



n 



E2 



"Ac 



d ^A(z) 
dz2 



+ u?A (z) 



(A20) 



^^ ^ic - rip t^R l^rPoL 

Now, we diff"erentiate Equation (IA20I) with respect to z twice and compare the resulting expression 
with Equation (IA19I ). obtaining 



d^A(z) 
dz4 



^PpFip+PcV^ ^^ 



+ OJ 




d^A(z) 
dz2 



A (z) = 0. 



(A21) 
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The general Equation (65) of IDymova & RudermanI (|2006h and our Equation (IA21I) are equivalent 



if a constant piecewise density is assumed in the former and Cowling's diffusion is omitted in the 
latter. Equation (IA21I) can be solved by taking a solution of the form exp (ik^z) and obtaining the 
subsequent forth-order polynomial for L. Two independent values of k^ are possible, namely k^i 
and k^2- Thus, the general solution of Equation (IA21I) is 

A (z) = Di exp (ik^iz) + £>2 exp (-ik^iz) + D^ exp {ik^2z) + D^ exp (-ik^jz) , {All) 

with Di, D2, 1)3, and D4 constants that are determined by the boundary conditions at z = +Lp/1. 
However, to keep this general analysis implies that the following expressions are complicated and 
require an additional mathematical effort. Instead, we choose a more restrictive way to simplify 
matters. 

For our next analysis. Equation (|A21I) is rewritten in a convenient form as 

b'-d'A(z) d^A(z) (^\,, ^ ,..., 

0)^ dz^ dz^ c^p 

where c^ and b^ are defined in Equations (fT2l) and (fT3]) . 

In the case of the fundamental mode, one could assume that, when the terms related to Cowl- 
ing's diffusion and resonant absorption are present, the characteristic scale for the variations of the 
eigenfunctions in the z-direction is only slightly modified with respect to the ideal case without 
transitional layer. Therefore, a reasonable approximation is to relate the forth-order derivative of 
A (z) in Equation (|A23I) with the second-order derivative as follows 

where the quantity 7C plays the role of the longitudinal wavenumber. We can approximate 7C by 
its expression in the ideal case, namely 

"K^ « — , (A25) 

with c? the ideal kink speed. Hence, Equation (IA23I) becomes 



-kp 



''''*« + ^^Afe).0, (A26) 



dz- ~2 



«u> .-, 



kr,/ 



which is formally identical to Equation (IA16I) . It is important recalling that the approximation 
of the forth-order z-derivative of A (z) may introduce some uncertainty in the solutions of Equa- 
tion (IA26I ) in comparison with the solutions of the full Equation (|A21I) . However, we expect a 
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minor discrepancy in the case of the fundamental mode because its characteristic scale in the z- 
direction should not be essentially modified when the terms related to Cowling's diff"usion and 
resonant absorption are taken into account in the equations. 



The solution of Equation (|A26I) is 



A (z) = El cos 



CO 



v^^F^.) 



E2 sin 



Ci) 



V^ 



Cki 



if 



\z\ < Lp/2 



p/^5 



(A27) 



with El and E2 constants. When the prominence thread is centered within the magnetic tube, 
i.e., L~ = Lg = \{L- Lp j, the solutions of Equation (IA26I) can be separated according to their 
symmetry about z = 0. Thus, even modes are described by Ei i^ 0, E2 = 0, and odd modes by 
El = 0, £2 ^ 0. 



A.4. Matching the solutions at z = +Lp/2 



In order to match the solution in the prominence thread (Equation (IA27I) ) with those in the 
evacuated regions (Equation (lAni) ). we impose the boundary conditions 



[[A]] = 0, 



dA 
d7 



0, 



at 



+LJ2 



p/^j 



(A28) 



corresponding to a contact discontinuity (JGoedbloed & Poedtsll2004l) . After applying these bound- 
ary conditions, the general dispersion relation (Equation (ITTI) ) is finally obtained. 
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